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We report the first observation of the radiative decay B° — > <j>K j using a data sample of 
772 x 10 6 BB pairs collected at the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. We observe a signal of 37 ± 8 events with a significance of 
5.4 standard deviations including systematic uncertainties. The measured branching fraction is 
B(B° — ¥ <f>K°^j) = (2.74±0.60±0.32) x 10 -6 , where the uncertainties are statistical and systematic, 
respectively. We also report the first measurements of time-dependent CP violation parameters: 
S 4,k%~, = +0-74±!J;^ (stat)to.2°(syst) and A^o^ = +0.35 ± 0.58(stat)lg;^(syst). Furthermore, we 
measure B(B + -> <f>K + ~f) = (2.48 ± 0.30 ± 0.24) x 10~ 6 , Acp = -0.03 ± 0.11 ± 0.08 and find that 
the signal is concentrated in the M^k mass region near threshold. 

PACS numbers: 14.40.Nd, 13.25.Hw, 11.30.Er 



Rare radiative B meson decays play an important role 
in the search for physics beyond the standard model 
(SM). These are flavor changing neutral current decays, 
forbidden at tree level in the SM but allowed through 
electroweak loop processes. The loop can be mediated by 
non-SM particles (for example, charged Higgs or SUSY 
particles) and therefore is sensitive to new physics (NP). 
Here we report the first observation of a new b — )■ s 
radiative penguin decay mode, B° — > (f>K°j, as well 
as measurements of its time-dependent CP asymme- 
try. This type of decay is sensitive to NP from right- 
handed currents [l[ and will be useful for precise time- 
dependent measurements at future high-luminosity flavor 
facilities 0-0]. 

The emitted photons in b — >• (b — > sj) decays are 
predominantly left-handed (right-handed) in the SM, and 
hence the time-dependent CP asymmetry is suppressed 
by the quark mass ratio (2m s /mb). The expected mixing- 
induced CP asymmetry parameter (S) is 0(3%) and 
the direct CP asymmetry parameter (A) is ~ 0.6% 
In several extensions of the SM, both photon helicitics 
can contribute to the decay. Therefore, any significantly 
larger CP asymmetry would be clear evidence for NP. 
In contrast to B° K*°(^ K s ir°)~f 0,0, another re- 
lated mode that is sensitive to NP, the time dependence 
of B° — > (fiKg'-f can be measured from the <f) ~^ K + K~ 



decay and docs not require a difficult measurement of the 
long lived K s decay inside the inner tracking volume or 
reconstruction of a low energy ir°. The B — > (f>Kj mode 
can be used to search for a possible contribution from 
kaonic resonances decaying to <pK . Furthermore, we can 
also probe the photon polarization using the angular dis- 
tributions of the final state hadrons 0] ■ 

Results on B — > 4>Kj decays have been reported by 
both Belle and BaBar collaborations based on 96 x 10 6 
BB @] and 228 x 10 6 BB pairs, respectively. Only up- 
per limits on B(B° — > (j>K 0r y) were given. Here we report 
the first observation of B° — > (fiK ^, the first measure- 
ments of time-dependent CP violation in this mode, as 
well as more precise measurements of B + — > cj>K + j [1 01 ] - 
The data set used consists of 772 x 10 6 BB pairs col- 
lected at the T(45) resonance with the Belle detector [ll| 
at the KEKB asymmetric-energy e + e~ (3.5 on 8.0 GeV) 



collider [12 1 



At KEKB, the T(4S') is produced with a Lorentz boost 
of j3j = 0.425 along the z axis, which is defined as op- 
posite to the e + beam direction. In the decay chain 
T(4S') — > B°B° —> / r oc.ftag, where one of the B mesons 
decays at time i rec to the signal mode / rec and the other 
decays at time £t ag to a final state /t ag that distinguishes 
between B° and B°, the decay rate has a time depen- 
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dence given by 



V(At) 



;-|At|/r B0 
4r R o 



i + q 



S sm(ArridAt) 



+ _4cos(Am rf Ai) j. (1) 

Here r^o is the neutral B lifetime, Am^ is the mass 
difference between the two neutral B mass eigenstates, 
At = t rcc — < t ag; and the 6- flavor charge g equals +1 (— 1) 
when the tagging B meson is a B° (B°). Since the B° 
and B° are approximately at rest in the T(4S I ) center- 
of-mass system (cms), At can be determined from Az, 
the displacement in z between the two decay vertices: 
At ~ Az/(Pyc). 

Signal candidates are reconstructed in the B + —> 
(f)K + -f and B° -» <f>K s j modes, with 4> -> K + K~~ and 
Kg — > ir + ir~ . Charged kaons arc identified by requiring a 
likelihood ratio Ck/tt [= Ck/{Ck + £■■*)] > 0.6, which is 
calculated using information from the aerogel Chercnkov, 
time-of- flight, and drift chamber detectors. This require- 
ment has an efficiency of 90% for kaons and an 8% pion 
fake rate. A less restrictive selection Ck/tt > 0.4 is ap- 
plied to the kaon candidates that are used to reconstruct 
the <j) meson. The invariant mass of the </> candidates is 
required to satisfy \M K + K - — < 10 MeV/c 2 , where 
denotes the <j> meson world-average mass [13(. The 
Kg selection criteria are the same as those described in 
Ref. [HI ; the invariant mass of the pion pairs should be 
in the range M„+„- G [482,514] McV/c 2 . The high en- 
ergy prompt photons must lie in the barrel region of the 
calorimeter (ECL), have a cms energy E^ ms e [1.4,3.4] 
GeV and a shower shape consistent with that of a photon. 
We also suppress the background photons from TT°(ri) 
77 using a likelihood C K o(£ v ) < 0.25, as described 



in Ref. [15 1 . 



We combine a <j) meson candidate, a charged or neu- 
tral kaon candidate, and a radiative photon to form a 
B meson. B candidates are identified using two kine- 
matic variables: the energy difference AE = £^ ms — 
and the beam-energy-constrained mass Mb c = 



pems 
beam 



(p™ ls /c) 2 , where £^oam is the beam en- 
ergy in the cms, and Eg lls and p c g ns arc the cms en- 
ergy and momentum, respectively, of the reconstructed 
B candidate. In the Mb c calculation, the photon mo- 
mentum is replaced by (E™* m — E^' K S ) to improve its 
resolution. The candidates that satisfy the requirements 
A/ bc > 5.2 GeV/c 2 and |AJ5| < 0.3 GeV are retained 
for further analysis. Using Monte Carlo (MC) simula- 
tions, we find nearly 12% (3%) of signal events in the 
charged (neutral) mode have more than one B candi- 
date. In case of multiple candidates, we choose the best 
candidate based on a series of selection criteria, which 
depend on a x 2 variable formed using the candidate's 
4> mass (and the Kg mass in the neutral mode) as well 
as the highest E^ nls and the highest Ck/tt in the charged 



mode. For events with multiple candidates, this selection 
method chooses the correct B candidate for the charged 
(neutral) mode 57% (69%) of the time. 

The dominant background comes from e + e~ — > qq 
(q = u, d, s, c) continuum events. We use two event-shape 
variables (a Fisher discriminant formed from modified 
Fox- Wolfram moments [l6| and the cosine of the angle be- 
tween the B flight direction and the beam axis, cos#b, in 
the cms frame) to distinguish spherically symmetric BB 
events from the jet-like continuum background. From 
these variables we form a likelihood ratio, denoted by 
lZ s / b . We require 7Z s / b > 0.65, which removes 91% of 



the continuum while retaining 76% of the signal. In ad- 
dition to the continuum, various BB background sources 
are also studied. In the B° —> 4>K S ^ mode, back- 
grounds from some b — > c decays such as D°7r°, D°r] 
and D~ p + , peak in the M\, c distribution. We remove 
the dominant peaking backgrounds by applying a veto 
to ipK^g combinations consistent within detector resolu- 
tion (±4cr) with the nominal D mass Some of the 
charmless backgrounds, where the B meson decays to 
<f>K*(892), 4>Ktt q and cj)Kr] also peak in M hc but shift 
towards lower AE. Another significant background is 
non-resonant (NR) B —> K + K~ Kj, which peaks in the 
AS- Mbc signal region; it is estimated using the </> mass 
sideband, M K+K - 6 [1.05,1.30] GeV/c 2 , in data. 
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FIG. 1: AE and M bc projections for B + — > <j>K + 'y (upper) 
and B a — > 4>KsJ (lower). The AE projections include the 
requirement M bc £ [5.27, 5.29] GeV/c 2 while the M bc projec- 
tions require AE € [—0.08, 0.05] GeV. The points with error 
bars are the data. The curves show the total fit function (solid 
red), total background function (long-dashed black), contin- 
uum component (dotted blue), the b — > c component (dashed- 
dotted green) and the non-resonant component as well as 
other charmless backgrounds (filled magenta histogram). 



The signal yield is obtained from an extended unbinned 
maximum-likelihood (UML) fit to the two-dimensional 
AS- Mb c distribution. We model the shape for the sig- 
nal component using the product of a Crystal Ball line 
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shape [13] for AE and a Gaussian for Mb c - The con- 
tinuum background is represented by the product of a 
first-order polynomial for AE and an ARGUS [l8| func- 
tion for Mb c . The b — > c background is described by the 
product of a second-order polynomial for AE and the 
sum of an ARGUS and a Gaussian function for M\, c . For 
the small charmless backgrounds (except for the NR com- 
ponent), we use the sum of two Gaussians for AE and 
a Gaussian for Mbc- The probability density function 
(PDF) is the product of these two functional forms fl9j |. 
In the final fit the continuum parameters are allowed 
to vary while all other background parameters are fixed 
to the values from MC simulations. The shapes of the 
b — > c and NR peaking background components are fixed 
to that of the signal. In the charged mode, the NR back- 
ground yield, (12.5 ± 6.7)% of the signal, is fixed from 
the </> mass sideband. Since the neutral mode is limited 
by statistics, we assume isospin symmetry and use the 
same NR fraction. The signal shapes are adjusted for 
small differences between MC simulations and data us- 
ing a 5° ^ i^*(892)°(-> A"+7T-)7 control sample, with 
M K+7r - G [820,970] MeV/c 2 . The fit yields a signal of 
144 ±17 5+^ 0A'+ 7 and 37 ± 8 B° -> ^A° 7 events. 
The projections of the fit results onto AE and M^ c are 
shown in Fig. [TJ The signal significance is defined as 
yf— 2 ln(£ /Anax), where £ max is the maximum likeli- 
hood for the nominal fit and Cq is the corresponding 
value with the signal yield fixed to zero. The additive 
sources of systematic uncertainty (described below) are 
included in the significance by varying each by its er- 
ror and taking the lowest significance. The signal in the 
charged mode has a significance of 9.6 o. whereas that 
for the neutral mode is 5.4 a. 
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FIG. 2: Background-subtracted and efficiency-corrected (j>K 
mass distributions for the charged (left) and neutral (right) 
modes. The points with error bars represent the data. The 
yield in each bin is obtained by the fitting procedure described 
in the text. A three-body phase-space model from MC simula- 
tion is shown by the filled circles (blue points) and normalized 
to the total data signal yield. 



To measure the M^k distribution, we repeat the fit 
in bins of <pK mass and the resulting signal yields are 
corrected for the detection efficiency. Nearly 72% of the 
signal events are concentrated in the low-mass region, 
M^k 6 [1.5,2.0] GeV/c 2 , as shown in Fig. H The MC 
efficiencies are reweighted according to this M^k depen- 



dence. These spectra are consistent with the expecta- 
tions from the pQCD model for non-resonant B — > 4>K^ 
decays [2(| ■ With the present statistics no clear evidence 
is found for the existence of a kaonic resonance decaying 
to 4>K. 

From the signal yield (A s j g ), we calculate the branch- 
ing fraction (£>) as A s ; g /(e x Ag-g- x £> sec ), where e is the 
weighted efficiency [(15.3 ± 0.1(stat))% for the charged 
mode and (10.0±0.1(stat))% for the neutral mode], N B -g 
is the number of BB pairs in the data sample, and jB sec 
is the product of daughter branching fractions [HJ . We 
obtain B(B+ ->• <f>K+i) = (2.48 ±0.30 ±0.24) x 10" 6 and 
B(B° ->■ 0A° 7 ) = (2.74 ± 0.60 ± 0.32) x 10~ 6 , where the 
uncertainties are statistical and systematic, respectively. 

We evaluate the systematic uncertainties on the signal 
yield by fitting the data with each fixed parameter varied 
by its ±1 tr error, and then taking the quadratic sum of 
all differences from the nominal value. The largest con- 
tribution of 8.0% arises from the NR yield. The other 
sources of systematic error are from charged track effi- 
ciency (~ 1.1% per track), photon detection efficiency 
(2.4%), particle identification (1.4%), number of pro- 
duced BB pairs (1.4%), 4> and Kg branching fractions 
(1.2%), K s reconstruction (4.6%), and the requirement 
on lZ s /b (0.3%). The statistical uncertainty on the MC 
efficiency after reweighting is 1.0% (1.2%) in the charged 
(neutral) mode. Furthermore, we assign a systematic er- 
ror of 0.2% (2.7%) for possible fit bias, which is obtained 
from ensemble tests with MC pseudo-experiments. The 
total systematic uncertainty on the branching fraction is 
9.5% (11.7%). 

For the CP asymmetry fit, we select events in the 
signal region defined as M^ c G [5.27,5.29] GeV/c 2 and 
AE G [—0.2, 0.1] GcV. Different selection criteria on TZ B /t 
arc used depending upon the flavor-tagging information. 
In addition, ECL endcap region photons are included in 
the analysis. We use a flavor tagging algorithm [2l| to 
obtain the 6-flavor charge q and a tagging quality factor 
r G [0, 1]. The value r = signifies no flavor discrimina- 
tion while r = 1 implies unambiguous flavor assignment. 
The data are divided into seven r intervals. The vertex 
position for the / rec decay is reconstructed using the two 
kaon tracks from the <f> meson and that of the /t ag decay 
is from well-reconstructed tracks that are not assigned 
t° /rec [22^ . The typical vertex reconstruction efficiency 
(z resolution) is 96% (115 /zm) for / rcc and 94% (104 /xm) 
for /tag- After all selection criteria are applied, we obtain 
75 (436) events in the signal region for the CP fit with a 
purity of 45% (37%) in the neutral (charged) mode. 

We determine S and A by performing an UML fit to 
the observed At distribution by maximizing the likeli- 
hood function C(S,A) = Yli Pii^i A] At,), where the 
product is over all events in the signal region. The like- 
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lihood Pi for each event is given by 
Pi = (1 - /oi) / hr ./Vr.i A/' ;/,'...: A/, - Ai')j d(Af) 

+ / iP i(A^), (2) 

where j runs over signal and all background components. 
Vj(At) is the corresponding PDF and Rj(At) is the At 
resolution function. The fraction of each component (/_,-) 
depends on the r region and is calculated for each event 
as a function of AE and M^ c . The signal PDF is given by 
a modified form of Eq. ffl) by fixing rgo and Amd to their 
world-average values [131 ] and incorporating the effect of 
incorrect flavor assignment. The distribution is then con- 
volved with a resolution function to take into account the 
finite vertex resolution. Since the NR component is ex- 
pected to have the same NP as the signal B — > 4>K^f, we 
treat this as signal for the time-dependent fit [23| . For the 
other BB components, we use the same functional forms 
as signal with an effective lifetime taken from MC and 
CP parameters fixed to zero. For the continuum back- 
ground, we use the functional form described in Ref. [22| ; 
the parameters are determined from a fit to the At distri- 
bution of events in the data sideband M^ c < 5.26 GeV/c 2 
and AE G [0.1,0.3] GcV. The term P Q i(Ai) is a broad 
Gaussian function that represents an outlier component 
with a small fraction f Q \. The PDFs and resolution func- 
tions are described in detail elsewhere pjj . 

We perform various consistency checks of the CP fit- 
ting technique. A lifetime fit to the B a — > A'*°(-> 
K + tt~)j, B + -> (f>K + j and B° <j)K%^ data sample 
yields 1.56 ± 0.03 ps, 1.70 ± 0.20 ps and 2.09 ± 0.45 ps, 
respectively. These are all consistent with the world- 
average values of the B lifetimes. The results of the 
CP asymmetry fit to the B° -)■ A"*°(^ (S = 

+0.02 ± 0.06, A = -0.06 ± 0.04) and B+ — » 0A+7 
(S = +0.25 ±0.33, A = +0.18 + 0.26) are consistent with 
zero. A fit to the sideband events in the B° — > <pKgj 
data sample gives an asymmetry consistent with zero 
(S = -1.77 ± 1.30, A = -0.04 ± 0.14). 



The results of the fit arc S = +0.74±?;^(stat)tg;^(syst) 



+0.35 ± 0.58(stat)Io;^(syst), where the uncer- 




FIG. 3: At distributions for q = +1 and q — —1 (left) and 
the raw asymmetry (right) for well-tagged events. The dashed 
curves in the At plot are the sum of backgrounds while the 
solid curves are the sum of signal and backgrounds. The solid 
curve in the asymmetry plot shows the result of the UML fit. 



The only free parameters in the CP fit are S and A. 



and„4- Ln«xncQ^+0.23/ 
tainties are obtained as described below. We define the 
raw asymmetry in each At bin by (N + — AL. ) / (N + + A_ ) , 
where N + (A_) is the number of events with q = +1 
(— 1). Figure [3] shows the At distributions and raw asym- 
metry for events with good tagging quality (r > 0.5, 48% 
of the total). 

We find that the error on S in the MINUIT minimiza- 
tion [24j is much smaller than the expectation from MC 



simulations and has a probability of only 0.6% [25J. This 
is due to low statistics and the presence of a single spe- 
cial event (with At = —3.64 ps and r = 0.96). A similar 
effect was found in our early time-dependent analyses of 
B° -> 7T+7T- [II]. Instead of the errors from MINUIT, we 
use the ±68% confidence intervals in the residual distri- 
butions of S and A, determined from toy MC simulations 
as the statistical uncertainties on the result. 

We evaluate the systematic uncertainties from the fol- 
lowing sources. A significant contribution is from the ver- 
tex reconstruction (0.08 on S, 0.04 on ^4). We refit the 
data with each fixed parameter varied by its error to eval- 
uate the uncertainties due to signal and background frac- 
tions (0.03, 0.07), resolution function (0.02, 0.03), AA- 
M bc shapes (0.01, 0.01), continuum At PDF (0.01, 0.02), 
flavor tagging (0.01, 0.01) and effects of tagside interfer- 
ence 27J (0.004, 0.030). The uncertainty from physics 
parameters (t^o, Amd), effective lifetime and CP asym- 
metry of the BB background, is (0.05, 0.03). We also 
include a possible fit bias due to low statistics and the 
proximity of the central value to the physical boundary 
MC simulations show that this bias de- 
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creases to 0.04 with twice the signal yield. Adding all 
these contributions in quadrature, we obtain a system- 
atic error of ILq'^I on S and ILq'iq on A. 

In summary, we report the first observation of a new 
radiative decay mode, B° — > (pK ^ using a data sam- 
ple of 772 x 10 6 BB pairs. The observed signal yield 
is 37 ± 8 with a significance of 5.4 a including system- 
atic uncertainties, and the measured branching fraction 
is B(B° ->■ (j)K°j) = (2.74 ±0.60 ±0.32) x 10~ 6 . We also 
measure B{B+ -» (j)K+j) = (2.48 ± 0.30 ± 0.24) x 10~ 6 
with a significance of 9.6 a. Furthermore, we measure the 
charge asymmetry A C p = [N(B~) - N(B+)]/[N(B-) + 
N(B+)} = -0.03 + 0.11 + 0.08, where N(B~) and N(B+) 
are the signal yields for B~ and B + decays, respec- 
tively. The signal events are mostly concentrated at 
low <j)K mass near threshold. The branching fractions 
and 4>K mass spectra are in agreement with the theo- 
retical prediction of Ref. We also report the first 
measurements of time-dependent CP violation param- 



eters in the neutral mode: S = +0.74 



+0.72+0.10 
-1.05-0.24 



and 

A = +0.35 ± 0.58i^;io- We have established that the 
mode B° — > ipKgj can be used at future high luminosity 
e+e ~ S 13 an d hadronic facilities Q to perform time- 
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dependent CP violation measurements and to carry out 
sensitive tests for NP. 
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